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Abstract 

New constructions in the theory of fields for multiple integrals are designed. 
Generalizations of the Legendre - Weyl - Caratheodory transforms and correspond- 
ing invariant integrals are introduced and explored. Connection and curvature of 

■ bundles induced by a field of extremals are calculated. 

(N 
H 

1 Introduction 

Let 91 be a domain on a smooth n -dimensional Riemannian manifold, and let p : £ — > 9T 
be a //-dimensional vector bundle over the base 9t. A fiber of the bundle over a point 
t e 9t, i.e. the full inverse image of the point t under the map p , is an v -dimensional 
linear subspace. Local coordinates on 91 are denoted by t = (t 1 , ...t n ) ; local coordinates 
on fibers are denoted by u ) . Usual convention on summation over two-fold 

occurring indices is used. Latin indices correspond to coordinates on the base and vary 
from 1 to n , Greek ones correspond to coordinates on fibers and vary from 1 to v . 
Multi-indices will be denoted by capital Latin and Greek indices respectively. Collection 
of all indices from 1 to n is denoted by / . The ordered exterior product of differential 
entering into the multi-index K is denoted by dt K := dt n A dtf 2 A ...dt lk (the letter K 
(N) \ shows that \K\ = k, where the absolute value of multi-index means the number of its 
indices). The symbol wedge for exterior product (the symbol A) will be omitted for the 

■ sake of brevity; it always will be implied while dealing with product of differentials. 
Consider the functional whose part related to the chart V C 9T , is 

Subsequent calculations will be produced in coordinates of the chart V . 

■ Let us denote by Ji(£) the bundle of 1-jets over £ and let 

a i 

f)T a 

<!?■■= w = g?(t,x) (2) 

be a section of Ji(£) . Such section can be considered as a slope field Q, i,e. as a 
distribution of n -dimensional planes in the space £ . We shall say that a manifold 
M = {x = x(t)} embedded in the slope field Q , if ^ = gf(t, x) . 
We distinguish three standpoints concerning the arguments in our functions: 

1. If t l ,x a , are taken as independent variables, as for instance in the function L , 
then the derivatives with respect to these variables are marked by attaching the respective 
variable as an index. 

2. By using a given slope field Q , the arguments ^ are replaced by functions 
g?(t,x) that depends on the f and x a only. The partial derivatives with respect to the 
arguments t l and x a are then denoted by ^ and . 
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3. The substitution x(t) and its derivatives referring to a given surface V changes 
functions which appeared in the second (or the first) standpoint, into functions of the t 
alone. Their derivation with respect to P is denoted by . 

The vanishing of the first variation is expressed by Euler's equations 

df dx a 
dP J 

Functions x(-) which satisfy this equation are called extremals. The Jacobi matrix 

of variables S relative to the arguments t is denoted by -ppf • The calligraphical V 
means the derivative with respect to the explicitly entering argument, while the direct 
D (say "Tjjpy- ) means the full derivative taking into account the dependance x{t) . The 

determinant of this matrix is denoted by | ^mpr I • The identity matrix is denoted by I 
(its dimension is implicitly defined by the formula in question). 

Nonnegativity of the second variation is the natural necessary condition of mini- 
mum |8j. The investigation of conditions of the nonnegativity begins with works of 
A.Clebsch JT3] who explored the Dirichlet functional 



b 2 T 



d 2 f 



dh a dh? 



+ 2- 



d 2 f dh a 



dP 



d 2 f 



dx a dx /S 



dt 1 . 



The idea of Clebsch was to reduce this functional to the integral from its main part, i.e. 
to the rearrange quadratic terms relative to first derivatives. The reduction was realized 
by adding under the integral sign a closed differential form having the type of divergency. 
It seems that Clebsch presumes that for multiple integral a direct analog of Legendre 
condition is valid: if the second variation is nonnegative then the quadratic terms relative 
to first derivatives that is defined on the space of (n x v) -matrices 



dx c 
~dP 



must be nonnegative for all values of t . 

But J.Hadamard [20] almost through half a year after the work of Clebsch shows that 
it is not correct. This quadratic form is nonnegative not for all matrices. The correct 
necessary condition (which is called Hadamard-Legendre condition) is the following. 



Theorem 1 (Hadamard) . Let the functional 



5 2 F 



v 



dP dti 



dt 1 



(3) 



be nonnegative for functions x(-) that meet boundary conditions 



x\ dv = 0. 

Then for all values of t £ V the quadratic form a l ^Jt)q"qj takes nonnegative values 



on in x v) 



-matrices having the form qf = ^ a r]i (that is for matrices of rank 1) 
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The assertion of the theorem can be reformulated as follows: the biquadratic form 
a ap( t )€ a € Pr 1i r lj is nonnegative for all values of t G V and £ eR u ,i] 6 (M n )* . 

Let us remark that a simple sufficient condition for optimality of small pieces of ex- 
tremals is the condition of convexity of / relative to variables Dx/Dt. But this assump- 
tion is much more strong than the necessary condition of Hadamard-Legendre. 

The work of Hadamard stimulates whole series of works that aim for decreasing the 
gap between necessary and sufficient conditions for optimality. One of the bright work of 
this cycle was that of Terpstra [30J where abstract algebraic questions suggested by this 
themes were considered. 

Namely, let 

be a quadratic form on {n x v) -matrices Wq^W ■ Consider the cone of rank 1 matrices, that 
is ones having the form qf = £ a (g) 77, . On this cone the form (j3j) turns into a biquadratic 
form 

4*(oreW: (5) 

which is defined on pairs of vectors £ G R", rj e (R n )* . Assume that the form is 
nonnegative on matrices of rank 1. Terpstra set up the question: is it possible to turn (j4]) 
into a positive form on the space of all matrices by adding terms 

where tensor r^g is skew-symmetric relative to both and a, f37 Similar additions 
to the main terms arise by adding a closed differential form to the integrand. These 
additional terms give zero on rank 1 matrices, and the biquadratic form (J5D do not changes. 
Terpstra shows that it is possible under the following condition 

Condition 1. The form (TJP admits decomposition into a sum of squares of bilinear forms 
from which as a minimum nv forms are linearly independent. 

The question on a decomposition of even forms into a sum of squares was investigated 
by Hilbert. The final result is due to Van-der-Waerden who proves that a decomposition 
is realizable if min{y, n) < 2 . If this minimum is more or equal to 3 then there exist 
indecomposable positively definite forms. The first explicit example of such forms was 
constructed by Terpstra [3U]- More simple example was suggested later by D.Serre [2~S] . 



This subject appears closely related with problems of existence of minima. The leading 
part in the proofs of existence plays the condition of lower semicontinuity of a functional 
in question. 

C.Morrey [23] and later J. Ball [9] show that necessary and sufficient condition for 
lower semicontinuity of second variation is the Hadamard-Legendre condition. And they 
explored an interesting generalization of the notion of convexity for multiple integral that 
is called polyconvexity. To define it, consider an integrand f(t, function of 

matrix 

dx a 



dt} 
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where t and x are fixed. One corresponds to each (n x v) -matrix the set of elements 
of its exterior powers, that is all (I x /) -minors of the matrix 1 < I < min(n, v) . One 
obtains a point of r -dimensional space denoted by r(g) . It is easy to calculate that 
r = ( n+ "j — 1 . The mapping r transfers M nv into an algebraic subset K, of the space 
W which is defined by Pliicker relations on minors of the matrix q . Thus the function / 
appears to be defined on the Pliicker cone K of the space W . The function / is called 
polyconvex if it admits a convex extension on the whole space W . The polyconvexity of 
/ is a sufficient (but not a necessary) condition of lower semicontinuity of the integral (see 
[23] . [T9] . [26] . [TO]). It is interesting that functions constructed by Terpstra and Serre 
gives examples of lower semicontinuous but not polyconvex functionals. The gap between 
the necessary and the sufficient conditions for optimality was essentially shorten by Van- 
Hove [32] who proved that natural strengthening of the Hadamard-Legendre condition 

® 2 f tatp 



eevw > ^l 2 M 2 (e) 



gives the locally sufficient condition for C 1 -minimum. The expression "locally sufficient" 
means that the domain of the integration is sufficiently small. 

The idea of the Van-Hove's proof is the following. First, we make coefficients to be 
frozen, i.e. we fixe arguments (t = to,x = Xq) in coefficients of the quadratic form of the 
integrand ([3]). It does not affects on the estimations since it can be taken as a domain of 
integration an arbitrarily small neighborhood of the point (t = to,x = Xq) . Second, we 
apply the Fourier transform to the functional obtained and use Parseval equality. 

This construction disclose the internal reason why the Hadamard-Legendre condition 
includes only rank 1 matrices. Namely, Fourier transform transfers the operation of differ- 
entiation into the operation of multiplication by the corresponding independent variable. 
So, if Fourier image of a function x a (t) is £ Q (7?) then the image of its derivative will 
be £, a r]i . As a result, the biquadratic form which stand at the left side of the formula (JB} 
appears under the integral sign of the Parseval equality, and the inequality (JBJ guarantee 
the positive definiteness of the functional in question. 

To prove optimality conditions on "large" parts of the manifold, one needs the theory 
of index of the functional, and the generalization of the notion of the conjugate point. It 
was the subject of big series of works (see, for example, Dennemeyer [IB] . Simons [27] . 
Smale [28], Swanson [29], Uhlenbeck [31]). We do not concern this theme here. 

The positivity of the second variation is insufficient to prove the sufficient conditions 
of the strong minimum. One needs the approach connected with the field theory and 
invariant Hilbert integral [18], [IT]. The first variant of this theory was suggested by 
C.Caratheodory [12J. Another (relatively more simple) variant with more strong demands 
on functions was constructed by H.Weyl [53]. We will shortly describe both approaches 
in a suitable form for subsequent presentation. 

In variational calculus for multiple integrals one needs as many principal functions of 
Hamilton (action-functions) as there are independent variables. Let we have n action- 
functions: S l (t,x), (i = 1, ...n) . The Weyl construction is based on invariant integral 
using divergence of a vector S l : 



Q = \"dt l ...dt i - 1 dS i dt i+1 ...dt n = —dt 1 . 

^ dP 

i=i 
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Let us recall that under the expression ^ we mean the full derivative of the function 
S 1 taking into account the dependence x = x(t) . For functions of matrices qf = dx a /df 
Weyl suggests to use the direct analog of Legendre transform in which the part of the 
scalar product plays the trace of the product of matrices. 

f(t,x,q) ^ f*(t,x,p*) = -/ + tr(p*g), (p% = J^. 

Here the value q as a function of p* has to give global maximum to the Weierstrass 
function — / + ti(p*q) . So one needs the condition of convexity of the function / on 
the space of matrices q with fixed t and x . The condition for the minimum of the 
Weierstrass function is an equality 

dS i df 



dx* 9(f£) 



The construction of Caratheodory is based on an invariant integral of the determinant 
from derivatives of functions 5** . 



n = det 



dS i 



dP 



dt 1 . 



Caratheodory considered the Legendre transform as an algebraic mapping of the space 
of quadruples. Let (f,(f>,q^,p z a ) 1 < i < n, 1 < a < z/ be a set, consisting of 2{nv + 1) 
elements, where / is the initial function, is the dual function, q® are the initial 
independent variables that correspond to "velocities" ( dx a /dt % ), and p l a are the dual 
variables that correspond to moments (df /d(^-)) . It is supposed that all these variables 
are bound together by the relation 



/ + = tr(pg). (7) 

We shall speak that a transformation is tangential if it transfers the canonical dif- 
ferential form df — pdq on the bundle of 1-jets into proportional to it canonical form 
df* - q*dp* . 

New variables of the standard Legendre transform are obtained from the old ones by 
simple permutation 

/* = 0, <P* = /> Q* =P, P* = Q- 

It is evident that this transformation is tangential, birational, involutory, and it pre- 
serves the relation (jTJ). 

Caratheodory has constructed a new transformation in the space of quadruples for 
the theory of invariant integrals having the type of determinant. The original text of 
Caratheodory is difficult, because of lack of motivations. Big series of many cumbersome 
formulas with numerous indices, which at the end bring to desirable results, give the 
impression of magic. So, we will give here the proof of the main theorem of Caratheodory 
in convenient for us matrix-form. 

Let us introduce an auxiliary square matrix A 

A=\\fl-pq\\. (8) 
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Introduce the notation f-r-j = 7 • The transpose of a matrix B will be denoted by 
Define the mapping Z of a quadruple (f,<fi,p,q) by the formulae 

f* = if; 0* = 70; ; (p*Y = iqA; (g*)* = A~ l p. (9) 

It follows 

p*q* = ApqA^ 1 ^. 

Using permutability of matrices A and pq , one has 

p*q* pq 

~F = 7 

Hence, taking into account the condition f + <f> = tr (pq) , one obtains 

r + 0* = tr(pV)- 

Theorem 2 (Caratheodory). . 

The transformation Z is tangential, birational, involutory, and it preserves the rela- 
tion (fjj. 

Proof. 

Rewrite the formula (E]) as pq = fl — A and express (q*) l q using ([9]). One has 
(q*Yq = A~ l pq = A~\fl - A) = A~ l f - I . Hence, I + [q*fq = /A" 1 . It follows 

£n 

det(I+ (q*) t q) = -r^ J = fF. (10) 

Calculate the differential of (11 01) using the following formula for differential of deter- 
minant 



We have 



d(det g) = (det g) ti(g dg). 

fdf* + f*df = ff* (tr( j[(dq*Yq + (q*ydq] ) ) = 

r (tr(^f^)+tT(A( q *Yd q )) = (n) 

/r(tr(^l)+tr(f )) = 
ftr(p*dq*) + f*ti(pdq). 
By combining the first and the last members in the chain ffTTl) . one has 

f(df* ~ tr(p*dq*) + f\df - tripdq)) = 0. (12) 

Now let us prove that the transformation Z is birational and involutary. Indeed. 
Repeat the transformation Z o Z . Introduce the auxiliary matrix A* = f*I — p*q* . 
Using ([8]) one obtains 

A* = 7/I - {iqA)\A- l py = 7 ||/I - {A- l pqAy\\. 
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The matrix pq commute with the matrix A , hence 



(/*) 



,n— 2 fn—2 



A* = 1 \\fl-(pq) 



* 
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det A* 



7™ det A 



1/7- 



(13) 



Using (jUJ), fTTH]) successively obtain 



p** = 7 *( g *A*)* = 7*7AA _1 p = 



= g- 



(14) 



□ 



Later on Th.De Donde [15], J.T.Lepage |21j . H.Borner [TT] a.o. connect approaches of 
H.Weyl and C.Caratheodory with the technique of differential forms. They tried to find 
the most general differential forms giving the field theory. 

We think that one has to find only invariant differential forms having direct geometrical 
meaning (such as forms of Weyl and Caratheodory). Action- functions define simultaneous 
parametrization of all the manifolds which define "the flow" of solutions generated a field 
of extremals. This is the reason why differential forms that characterize "the flow" in the 
bundle £ have to be invariant relative to choice of coordinates for a mapping £ — > S . 
Besides, it is not imperative that unknown manifolds admit one-to-one projection onto 
the space of variables (t 1 , ...t n ) . This fact counts in favour of invariance. 

The main goal of this work is to find and to explore new constructions of invariant 
integrals and the corresponding tangential transformations. It will be designed a series of 
transformations using tangent planes of f(Dx/Dt) in spaces of different exterior powers 
of matrices of jets. Note that Caratheodory uses only one (the highest) exterior power - 
the determinant. 

To construct a field theory one has to add to an integrand a closed differential form, 
that does not changes the value of the integral but turns the integrand into nonnegative 
function (analog of the Weierstrass function). As in the theory of characteristic classes [6], 
to design differential forms it will be reasonable to use invariant polynomials on the matrix 
algebra. Such polynomial were the trace and the determinant, i.e. forms suggested by 
Weyl and Caratheodory respectively. We consider remaining generators of the ring of 
invariant symmetrical polynomials from roots of matrices ||^|| • 




Let C be a square (n x n) -matrix. 
Definition 1. . 

The minor of a matrix C is called the principal one if the indices of its rows coincide 
with indices of its columns. 

Lemma 1. . 

The sum of all principal minors having the order k for a matrix C is equal to the 
coefficient of the characteristic polynomial -P(A) standing before the term \ n ~ k . 




2 k-Lagrangian Manifolds 
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Proof. 

It is sufficient only to see on the matrix 

C= ||4-A5;.|| 

and to note that the coefficient before \ n ~ k is obtained by elimination rows and columns 
of any k diagonal elements (that leads to principal minors), and subsequent summation 
over all such elements. 
□ 

Suppose that n action-functions S l (t,x), (i = l,...n) are given. To construct an 
analog of a Hilbert integral, consider the following closed differential form 

& k = (-l) r{K) dS il ...dS ik dt I \ K , (15) 

\K\=k 

where K = (ii,...ik). By r(K) we denote the number of permutations needed to put 
members of the exterior product of differential forms dt n ...dP k dt I ^ K in ascending order. 
It is easy to calculate that 

r i K ) = 2^ l » 2 — ' 

is&K 

The differential form & k corresponds to the coefficient of \ n ~ k of the characteristic 
polynomial -Pc(A) of the matrix 

OS 1 dS i dx a 
OP dx a dp 

Indeed, each of the summands of the matrix & k has the following Jacobian that 
corresponds to the principal minor of the matrix C as a coefficient 

D[S K ] 

For k = 1 we have Weyl theory, for k = n Caratheodory theory. 

To represent the differential form & k in terms of the integrand / we need a small 
regression. 

Consider (n x n) -matrix C = $ + \P , where $ = diag(0 fe ) is a diagonal matrix, 
\I/ = ||a l 6j|| is a rank 1 matrix being the tensor product of a contravariant vector a 1 
by a covariant vector bj . We will consider only the case (needed in what follows) when 
(p k , a 1 are differential forms of the first order, and bj are elements of the basic field. Due 
to the noncommutativity of the exterior product of differential forms, we will conceive 
determinant consisted in that elements in a special, not usual, sense. While expanding the 
determinant, we will order factors of each terms in accordance with numbers of its rows. 
Taking such rule, a determinant with two equal column (but not rows!) will be equal to 
zero. Hence, it is possible without changing of the determinant to add its columns (but 
one cannot add rows!). 



C 



dS l 
~d¥ 
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Lemma 2. The determinant of the matrix C is equal 

n n 

det C = Yl 0' + Yl ^■0 1 -0 J ~ 1 a i J+1 -0 n - 

i=l j=l 

Proof. 

Since among the components of the vector bj there exists at least one distinct from 
zero (otherwise det C = ), we can set without loss of generality that b\ ^ . Preserving 
the first column of the matrix C we subtract from each column with the number i the 
first one being multiplied by bi/b\ . The expansion of the matrix obtained relative to the 
first row gives the needed formulae. 

□ 

To find the analog of the Poincare-Cartan form, we rewrite integrand of (pQ) as follows 



_L_ £ ( _ ir <*) J2det(diagfdt)dt^ = £ ((-1)™ ]J fdt) dt^ . 

\k)J \K\=k i£K \k)J \K\=k \ idK J 

Then we add to each factor, standing under the sign of product, the canonical form 
for the distribution (j2j): u a := dx a — Y^=i9f^ being multiplied by the corresponding 
momentum p l a := f q ? . The lemma 2 gives the reason to consider as a Poincare-Cartan 
form the following expression 



where 



A 



\K\=k 



(16) 



-\y{K) 



rn 

i£K 



f dtl + lT,P« (dx a ~gT(t,x)dt i ) 



o=l 



(17) 



Here ]~[ stands for exterior product of differential forms, and into the function / we 
substitute the distribution ([2]). In parentheses of the formulae (ITT)) only the summand 
gfdt 1 takes part, because other summands of the type g^dP give zero for j ^ i, since 
the term dp meet twice in the exterior product. 

The choice of differential form (IT6l) can be substantiated as follows. Each summand 
of the differential form dS K dt I \ K is a simple multivector of an order k . We can 
write it as an exterior product of one- dimensional forms. Using the invariance, we select 
in each factor, entering in dS K , independent summands and take it as dt % . The rest can 
be justify as follows. The differential form A is reduced to fdt on the integral manifold 
of the distribution (|2]). Besides, the derivatives of integrands of both functionals relative 
to qf coincide. These facts are provide by the choice of numerical coefficients (?) and 
n/k . It is necessary to speak about minimum of the Weierstrass function that will be 
built using the differential form ( TT6)) (see below the section "Weierstrass function"). 

Note 1. The exchange of the determinant by this product is equivalent, in essence, to the 
restriction of a symmetrical multilinear Ad -invariant form, defined on a Lie algebra, to 
the Cartan subalgebra. This restriction uniquely defines this form JT$. 
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Suppose that the distribution (F2J) is integrable. Then the manifold 971 C , 

dim 971 = v + n defined by equations ([2]) is fibred by n -dimensional fibers. Since 9JT 

(due to its definition) has one-to-one projection onto the space £ , the foliation on 9JT 
induces a foliation of £ by n -dimensional fibers $. 

Definition 2. . 

A manifold 9Jt C , dim9Jt = z/ + n is called k -Lagrangian if the restriction of 
the differential form A to 9H is the closed form. 

The differential form A lead to the natural generalization of the Legenre transform. 
It corresponds to describing of the function / from matrix variable as an envelope of 
the family of all tangent planes to the surface / : ~R nu — > R considered as a function of 
exterior k -power of its argument j§ . 

Theorem 3. Let the distribution (Tj|) be integrable, and the manifold M be k -Lagrangian. 
Then the fibers $ C £ are solutions to the Euler equations 

± fa) + U = 0. (18) 

Proof. 

Calculate the differential of the form 

A* = -iLfdt + Y^^ s dtix --- dtis - X i dxa - gZdt is )dt i ° + \..d?*]dt I \ K . 

\k) U-lJ a=l s=l 

Summands which contain more than one form u a are omitted in this formula be- 
cause after differentiation and subsequent substitution the distribution ([5} they give zero. 
While differentiating, we consider only summands of the type dx x dt I . The result will be 
presented as four groups. 

0i = jn^Q^xdx^dt 1 , 
02 = jjjf q « l ^xdx x dt I , 

03 = Eti(i)/c(-S) d*W 
04 = EjU j^^f q *dt i °dt il ...dt i °- 1 dx x dt i °+ 1 ...dt ik dt I \ K . 

The sums 0i and 02 are generated from differentiation with respect to x x of the 
summand -p^rfdt 1 . The sum 03 is generated from differentiation with respect to x x of 

[k) 

the factor u as . In the sum 04 stand the full derivatives with respect to t ls . From u a is 
retained only dx x . To order differentials in 04 , it is necessary to permute dt ts and dx x , 
so the sum 4 changes its sign. 

Further we have to sum over K the expression obtained. The expression <pi as 
well as J2k 02 consist of equal summands, and its number is equal to the number of 
groups K . Hence, the coefficients m are cancelled. Each summand in the expression 
Y2i 03 as well as in Y^,i 04 is me t as many times as there are groups that it contain. The 
number of such groups equal ClZi) ■ Hence, the coefficients are cancelled. After 

this reduction 2 + Y2k 03 = . The remaining summands 0i + 4 give 
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□ 

Note 2. For n = 1 the set of solutions to Euler equation is finite-dimensional. For n > 1 
there are functional freedom to choose solutions. These solutions could be combined in the 
"Lagrangian manifolds" variously. This is the reason why there are many different kinds 
of fields of extremals for the case of multiple integral. 



3 Weierstrass function 

Having invariant integral of the Hilbert type corresponding to the form ([16]) one can build 
the analog of the Weierstrass function for the distribution of the slope field g?(t,x) [I]. 
To have nonnegative Weierstrass function for all values of derivatives it is necessary that: 
first, the value of the integrand corresponding to Hilbert integral on manifolds imbedded 
into the slope field g^(t,x) was equal to the value of /, and second, its derivatives with 
respect to Dx/ Dt were equal to that of the function / . Canonical forms u a vanish on 
manifolds imbedded into the slope field in question. Hence, calculating differential one 
can ignore summands containing products of form to" = — gf(t, x) . Let us denote by 

/ the value of the function / after substitution of the slope field 

=/(«.*.£) - E g^<ie t ||/ I + ^K)||. (19) 

It is evident that the coefficient jnjp^i standing before the sum provides the condition 

f = f . Coefficients standing inside of determinants provide coincidence of derivatives. 
Indeed, let us expand the determinant in ffl9|) . Summands containing = ^jj- — g?(t, x) 
in the first degree arise if we take only one factor standing on diagonal. The derivative 
stands on the m place of diagonal, and its coefficient is |g™ . It does not depend on 
K . This coefficient meets as many times as there are multiindices K containing index 
m , that is (^-l) times. The total coefficient will be equal to (2-1) I (fc) = ~ • This gives 
the coincidence of derivatives of the function / and of the subtrahend in the formula (fl9l) . 
Hence, the following relations will be valid 

88 

£(t,x,g,g) = 0; (t,x,g,g) = 0. (20) 

d{%r) 

Definition 3. . A slope field gf(t,x) is called geodesic for the differential form &k in a 
domain V C £ if the minimum of the Weierstrass function fW\) is reached at each point 
(t,x)EV for d f t =g. 

Note 3. Let us recall that the Legendre transform in case of simple integrals and the Weyl 
transform in case of multiple integrals describe a function f(-,-,Dx/Dt) with the help of 
support planes to the graph f : M. nu — > R . By contrast, the transforms , defined below 
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in section 6, will describe f with the help of multilinear (relative to variables dx a /dt l ) 
support manifolds to the graph f : M. nu — > R . These support manifolds can be regarded as 
planes in the space of multivectors of dimension k . 

Theorem 4. . 

If a manifold x(-) is imbedded into the geodesic field <?f (t, x) in the domain V then the 
functional in question reaches on this manifold the minimal value relative to any manifold 
with the same boundary lying in the domain V . 

In addition, the local minimum relative to dx/dt gives the sufficient condition for the 
weak local minimum while the global minimum gives the sufficient condition for the strong 
minimum in V . 

Proof. 

The proof follows by the standard way from the invariance of the Hilbert integral and 
from the positivity of the Weierstrass function in the given domain. 
□ 

Note that in this theorem we do not suppose that the slope field gf (t, x) is integrable. 
Nevertheless, the imbedding into a geodesic field gives minimum for one individual man- 
ifold. If the field is integrable then there exists a manifold with the slope defined by 
the field which passes through each point of the domain V . Due to the theorem 2 it 
gives the minimum, and, a fortiori, it is the extremal. Hence, in the integrable case, the 
manifold x(-) is imbedded into a field of extremals. Through each point of the domain 
V passes one and only one n -dimensional extremal giving minimum to the functional. 
The n -parametric family of //-dimensional level surfaces of action-functions S l = Const 
transversally cross these extremals. (Transversality conditions will be found below in sec- 
tion 5.) The German classical literature on variational calculus awards to the described 
geometrical object the name "Perfect or complete picture" (eine vollstandige Figure). 



4 Condition of solvability 

Let us find conditions for minimum of function £ at a point g . The first derivative at a 
point g must be zero. We differentiate £ by dx x /dt m . This argument enters only into 
the m -column of each determinant. The derivative will be equal to the same determinant 
in which the m -column is changed by 



By expanding the determinant relative to this column one obtains 

^-T^i-r E E(*))^) = ' ( 21 ) 

\k)J \K\=k,K^li&K 

Denote by C(K) the matrix corresponding to the principal minor with the index K 
in formula ( TT9|) . Denote by adj'(C(i^)) the adjunct (the cofactor) of the element (z, I) 
of the matrix C{K) . Denote by adj'™(C(i^)) the adjunct of the minor standing in the 
columns (l,m) and in the rows (i,j) of the matrix C(K) . 



12 



The second derivative of the function 8 at a point g is obtained by the exchange two 
columns by coefficients of the corresponding dx/dt. The expansion of these determinant 
relative to pair of columns (l,m) gives the following condition of minimum: The quadratic 
form with the coefficients 

gSL - £ £ % (^(c(K))) m - m w 

Hm \k) J \K\=k,K3(lm) (i,j)<=K 

on the space (nv x nv) -matrix must be nonnegative. Let us rearrange coefficients of 
quadratic form fl22|) which adds to the first summand 



d 2 f 



dq^dqifn 



(23) 



If the minor adj-™C(i^) is not the principal one ((/,m) 7^ (hj)) then, after the 



substitution q = g, it will have zero row or zero column. If (Z,m) = then it turns 

into the diagonal matrix with the diagonal elements / . So, we have adj*j C(K) = f k ~ 2 . 
Hence, all summands of coefficients in question appear the same. The number of these 
summands is equal to the number of minors K containing the pair of indices I, m . So 
the total coefficient is equal to (tlo) ■ Hence, we have 

a£^-|^ ( «-^ >0 - (24) 

The expression (p l x p™ — P^P™) defines the skew-symmetric form vanishing on matri- 
ces of rank 1. Hence, its addition do not violates Hadamard-Legendre condition. The 
formula (I24p gives precisely such skew-symmetrical forms which has to be added to the 
integrand to obtain the invariant integral of one or other degree k . 



5 Transversality condition 

Consider moving boundary problems of minimization of the functional ([T]). For example 
we can put that the boundary dV of unknown solution x(-) belongs to a fixed manifold 
X . It is evident that the solution to this problem meets the necessary conditions for 
optimality for the problem with the fixed boundary. Find additional necessary conditions 
caused by the possibility to vary the boundary. 

Suppose that each point of the manifold x(-) moves in space (t,x) along trajectories 
of a vector field T l (t, x),X a (t, x) that tangent to the manifold X . Denote by 9 the time 
of translation. The solution of the system 



(25) 



i = T(t,x) 
x = X(t, x) 

with initial conditions (to, x ) G x(-) will be denoted by T(9\ t Q , x ), X(9; t , x ) . By h(t) 
we denote the derivative of x with respect to parameter 9 for given t . 

dX 

h(t) = —(0;t,x(t)). 
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After substitution T(9),x(9) we obtain the function ^{9). Find 



d_ 

d9 



J-(O) 



dh a 

f x a h a + f dx<* ~ . 



dt 1 + / fT u dS. 



(26) 



Here T u is the projection of the vector T on u (on the normal to the boundary dV), 
and dS is the element of the volume of the boundary. Integrate by parts the summand 
including derivatives of h . Under the integral on V we get the left hand side of the Euler 
equation on the manifold x(-) that gives zero. Under the integral on dV it is added the 

summand /W; h a dt 1 ^ . Let us express h a at points dV through components of 

at 1 

the vector field (T,X) . 

Differentiation of the identity x(9) = x(T(9),9) gives 



dx a dx a dP 
~~39 + ¥¥ 



or 



Substitute fl2 
obtain 



h c 



into (|26l) and take into account that T u dS = T*(— 1 



,* — l 



dt 1 ^ 



(27) 
We 



d 

19 



J-(O) 



fdx^_X a + I /#!■ — f an" 
at 1 \ ■' at 1 



-1 



Theorem 5. T/ie necessary condition for optimality for moving boundary problem is 

( - ■ - dx a 

f dx a X a + I f8 l A — f dx a 



(28) 



c)V 



, , T ] = 

at 1 \ ' at ! ct J 

for any vector field T(t,x),X(t,x) which tangent to the manifold X . 
Proof. 

Let the theorem be violated for a vector T(t , xo), X(to, Xq) which tangent to X at 
a point (ioj^o) • Extend this vector to a smooth field which tangent to X and which is 
nonzero only in the sufficiently small neighborhood of the point (to, Xo) ■ For the variation 
corresponding to the shift along this vector field we get 

d 



Te m < o, 

that contradicts to minimality of x . 
□ 

The transversality conditions give n equations on the vector (T, X) . The number 
of conditions is equal to the number of action-functions. Clarify the meaning of these 
conditions in case k = n (the case of Caratheodory). We have 



fdt 1 



det 



dS i dS i dx a 
+ 



dP dx a dP 



dt 1 . 



Denote the matrix standing under the sign of determinant by F . Denote by a r s the 
matrix appearing in the transversality condition 
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ffi r _ fX^(F~ 1 Y 
J s J Z^V )i g xa dfs 



i=l 



The transversality condition is 



a r s T s + p r a X a = 0. 



In our case 



i=l 



- 1)r dS i dx a 
" dx a dt s 



n (99* 



1=1 



Divide by / and multiply from the left by the matrix F™ : 



dS m dx c 
dx a dt s 



a cm 

T s + - — X a = 0. 

dx a 



Substitute the value of F 



f) cm f) cm 

170 _|_ 170 -\ra _ g 

dt s dx a 

Consequently, the vector (T, X) is lying on the intersection of manifolds S m = Const. 



6 Generalization of Legendre- Weyl-Caratheodory trans- 
forms 



We fix value (t, x) as for the classical Legendre transform. These variables will not 
appears in subsequent formulas, and by f(q) we always understand f(t,x,q). Recall 
that 



dx a 

are the main arguments of the function / , and 

Of 



Pa 



are the corresponding moments. Denote by A k R exterior k -power of the matrix R . 
It is the matrix consisted from minors of order k of the matrix R. The function that 
corresponds to the invariant integral ffT5l) will be written in the form 



A 



Of"- 1 



trA A 



n 



(29) 



In view of lemma 1 it is natural to relate the expression A with the function 



6 = tr A A 



dS i <A dS i dx a 



a=l 



dx a dP 



(30) 
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The conjugate variables (with the variable qf ) which will play the part of new inde- 
pendent variable in our analog of Legenre transform must be tied with dS l /dx a . Keeping 
this in mind we introduce the matrix 



(lT m (t,x,q) 



dA 

0(qZ) 



(31) 



The differentiation of 



gives 



eans tl 

the auxiliary matrix 



n 



k 



n , 



\k)J K3(m,l) 

The subscript K means that entries of the matrix have indices belonging to K . Using 



[A 



E ( ad j< 



0/ 



K3(m,l) 



n 



f 5 ) + 



n 
k k 



we write in the form 



Denote the canonical function corresponding to / by 

H = f*(q*)= 1 (A-f). 



(32) 



(33) 



(34) 



The definition of this function will be given later in the formula (38). 

Let us express the differential form Q related to the function A in canonical coordi- 
nate. We carry out of brackets the function f*(q*) . Each simple multi- vector entering in 
the sum can be represented as a product of k differential forms of the first order. Select 
in each form the summand dt l , (i G K) . The number of these simple multivectors is (™) 
and we cancel Q by this factor. The coefficient of qf has to be equal to f q ? . Conse- 
quently, coefficients of dx a in factors of each monomial include (£Zi) " ~ the number of 
minors K containing the index i . In view of our cancellation by (?) this factor is | . 
As a result, under the integral appears multilinear function (relative to variables qf ). Its 
derivative with respect to gives (q*) 7 ? ■ Consequently, this function coincides with the 
function A . We arrive to the following differential form 



\K\=k \ieK ot=l 

Normalize variables by putting 



kH 



(35) 



q: 



kH 



The formula ( 135]) takes the form 



\K\=k \ieK a=l J 



(36) 



(37) 
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The condition of minimum of the Weierstrass function (j2ip can be written as the 
following equation p l a = dA/dqf , or q* = p t A' 1 . In view of ff24|) the equation (l3Tj) can 
be resolved relative to q . It gives the dependence q = 4>{q*) . 

Define the transform Zf. following the Caratheodory"s approach. Consider quadruples 
{/; 4>i qfiPa} > which satisfy the condition 



/ + 



A. 



Introduce the auxiliary matrix A by using ( 1321) . In view of ( 1241) the matrix A on 
extremals has to be nonsingular. The formula ( 1331) gives 



dA 
<9gm 

Denote by 7 the scalar coefficient 



7 



/ 



k-2 



detA 



Variables that are images of quadruples under the transform Zf. will be marked by 
superscript *. Define Zy. by formulas 



f = 70! 0* = 7/; (p*) = lAqS (q*) = V l A-\ 
The transformed function will be f*(q*) ■ 

Theorem 6. The transform : f(q) H- f*(q*) is birational and involutory. 
Proof. 

The fact that Zf. is birational is evident. 
The proof of involutority. 



(3t 



* * 

q p 



7 (p t A- 1 )(Ag i ) 

Since the transposition does not changes determinants we have 



A* = 7A. 

One obtains by multiplication of the formula / + 

/* + = A*. 



(39) 



A by 7 



The auxiliary matrix A 1 was defined through differentiation of the function ( 129]) 
Therefore the matrix (A^ 1 )* will be defined through differentiation of the matrix 



trA A 



n 



(r)^ + 7>W)j- 



In view of the formulas ( 13 8 p . (I39p . this matrix equals 



7 



(DP- 1 



trA A 



Tl 

f8] + - V \q«. 
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Hence A* = jA t , and 



« _ IT)*" 2 _ i k -\f) k - 2 _ i 

7 detA* 7 fe detA 7' 
The formulas obtained allow us to find the second iteration of the transpose Zf. : 

/** = 7 *(0*) = / ; 0** = 7 */* = 0; 

p ** = 7 *^*( g *)< = 7 * 7 A*(A- 1 )*p = p; (40) 
g « = (p*)*( i l*)- 1 = 7(^*)^(^- 1 )' = ?- 
The theorem is proved. 
□ 

In full agreement with the fact that Zu is involutory and tangential transform we have 
that the form Q is obtained from the differential form corresponding to the function (129]) 
by exchange all the variables by its dual (in the sense of the transform Zk ) . 

Note 4. There is an intimacy between the transforms Zk and the condition of poly con- 
vexity of function f . Indeed, if the function f admits a convex extension to the space 
V k (Dx/Dt) then it can be described as an envelope of support planes in the space of 
multivectors (see the note 3 in the section 3). Hence, after the addition of corresponding 
expression of the type (24\ ), we can apply to f the transform Zk . 



7 Formulas for action-functions 

Let us compare two expression of the integrand of invariant integral. The first is in terms 
of action-function 



6 



E det 

\K\=k 



dS i >A dS i dx c 



dP 



a=l 



dx a dti 



(41) 



The second is in terms of initial integrand / being written as integrand of differential 
form Q 




|if|=fc \ieK a=l 

We expand each of determinants of the formula (jUJ) in the sum of its columns and 
gather terms with the determinants of the same order. To write the corresponding for- 
mulas we put that J C K is a multiindex, and H C {l,2,...u} is a multiindex of the 
same order as the order of J . 



6 




V[S 



V[t K \ 



min(fc,^) 



E 



8=1 |J| = |S|=s 



V[S R \ J } 



V[S J ] 



V\x s ] 



V\ar] 



V[t" 



(43) 



Multiply brackets under the sign of the product in the formula (1361) . taking into account 

is equal to the coefficient arising by expressing of dx~ through 



that the Jacobian 



dt . Comparing 



V[tJ] 

with 



we obtain the series of formulas 
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E 

\K\=k 



E 



V[S 



A'l 



V[t K \ i x c 



V[S K ] 



H. 



Finally, to write the general formula we make denotation more concrete. We fix the 
sets J = {ji, ...j s } and H = {a%, ...a s } . Then 



E 

AdJ 

Expressions 



V[S 



A'l 



ff^det||Qi| 



AdJ 



i/ 1 " 8 det 



g[g 



are the Pliicker coordinates of k -multivectors composed from 



V[t K \ J x s ] 

gradients of the action-functions S l (relative to both dependent and independent argu- 
ments (x a ,t 1 )). Summarize: The canonical variables H and Q are sums of the Pliicker 
coordinates of gradients of the action-functions in the bundle £ . It is the usual gradient- 
vector in the Weyl construction, and it is only one multivector of maximum dimension n 
— the determinant — in the construction of Caratheodory. 

Now we find the corresponding expression using initial function / . Recall that 



A 



\K\=k \kJ J 



k-1 



det 



/3 + ?£&(#-#)) 



0=1 



(44) 



With the determinant in the formula (jUj) we shall carry out the same operation as 
with the function (I43p . We expand each determinant into the sum of its columns and 
gather determinants of the same order together. 

We calculate coefficients of the corresponding minors. Coefficient of / under the sign 
of the sum over K is equal (?) , and this sum include the same number of identical 
summands. So, the total coefficient equals to 1. Coefficient of f°p l a ^- under the sign 
of sum over K is equal to (?l!) and the sum over K includes the same number of 
identical summands, since it is just the number of fc-sets including the index i . So, the 
total coefficient equals to 1 too. We have CIZ2) summands including the pair of indices 

, and the coefficient of f~ l ^w3] * s equal to ^1 • The coefficient of the general 

- s+ lV\x E 



2?[f*>3] 

The minor composed from df/dp l a , where a G S, i £ J will be denoted by 



member /" 



We have the formula 



A 




-s+l 



(n) 
(ri) 



v[ P i] 



V[x E ] 



(45) 



If we equate coefficients of j^pj in formulas (Jilj) and f|45|) we obtain equations for 
action-functions in terms of / . 
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8 The canonical equations 



It is natural to write the analog of the Jacobi equation in the canonical coordinates, i.e. 
in terms of Pliicker cordinates of multivectors composed from gradients of the action- 
functions. To obtain these equations it is sufficient to write the conditions of closeness of 
the differential form Q : 

V = HJ2 {-l) r(K) [H(dt l + J2Qidx a )dt I \ K 

\K\=k \ieK o=l 

We have 

dn = z K (Bdt m + §^dx^(E K U(dt l + Q l a dx a )) 

+HJ2k (§*dx x + U(df + Q l a dx a ) = 0. (46) 

All the summands of this product are ordered relative to the order of dt % but in 
factors that do not contain the index i , i.e. i K , we omit all summands adding to 
dt % . Besides, in factors dt l + Q l dx 11 we omit that functions Q 1 ^ the derivatives of which 
was already calculated. We equate coefficients of the form dfl for independent products 
of differential. The coefficient of dx^dt 1 equals 



dH _ dH dQ 



III 



^ \dx^- Q ™d^- H m^ ]=0 - (47) 

KBrri x 

The sign "minus" for two last members in the formula ( H71) is stipulated by necessity to 
exchange the order of differentials dt m and dx 11 . 

All summands standing under the sign of the sum are identical, hence each of it is zero. 
To write equations in a more compact form, it is convenient to introduce an operator 



M dx» ^ M dt m ' 
In terms of the operator L the equation (H7|) takes the form 



H^- = LnH. 

dt m M 



Note 5. The abbreviated writing L a is not an occasional one. With the heuristic point 
of view, it is convenient to imagine Q l a as if it were the differential of t l with respect to 
x a do not calling attention to dimensions. Then the operator 

T ~ 9 -Q™ 9 



" 3x a ^ a dt m 

would be the full derivative with respect to x a taking into account the imaginary depen- 
dence t(x) . 

n\(q*y] m 

Rewrite the equation in an another form substituting in it Q™ = — kH M . We get 

n d[(g*)X _ n l( * V]m dH dH _ dH 
k dt m k m ' ^ dt m dx» Vm dt m 
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The second members in both sides of the last equation mutually annihilate. The equa- 
tion P8|) takes the form 



kdt m dx» ^ ' 

The definition f )34p of the function H and the involutivity of the trasform give 

d[( q *)t]m q ™ Qfm- ^ 

Combining ( )49l) and ( 1501) . we get the system 

n d[{q*)X _ dH 



dt m „ Qxi 1 



iT (51) 



Note 6. T/ie Euler equations were obtained by equating to zero the coefficient of dx x dt I 
of the differential form dA . The formula ( f7P| ) was also obtained by equating to zero the 
same coefficient of the same form but being written in canonical variables (in terms of 
Pliicker coordinates of gradients of the action-functions). Hence, it is natural to consider 
the system of equations / T51]) as a canonical form of the Euler equation. As in the classical 
case, the equation follows from very definition of the functions H and q* . The 
shape of the canonical system is the same for any k but the Hamiltonian as well as the 
function q* depends on the choice of the invariant integral. 



9 Necessary and sufficient conditions 
for closeness of the form fl 



Calculate the coefficient of dx ,1 dx x dt I \ l at the left hand side of the formula fj46|) . 

The summand ^^dt m gives two members. For one of these members we choose the 
differential dx^ 1 standing on the m-spot, and differential dx x on the 1-spot. To reorder 
we need / transpositions. For another one we choose the differential dx^ standing on the 
1-spot, and differential dx x on the m-spot. To reorder we need (I — 1) transpositions. 
These two members take the form 

fjlT 

(-l)^(^M-^)- 

The summand J^cfa^ gives two members too. By the same reasoning we verify that 
its sum is 



" ' \dx x ^» dx^ x 



The summand -g-^-dx from the second part of the formula (H6l) gives the following 
pair of members 
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■l) l H 



dQl_dQ[ 
dx x dx^ 



Finally, the summand -^dt m gives the fourth pair of members 



'_i v h I ;n£An m — — _ ^ n m 

For the fixed / the expression standing under the sign of the sum by multiindex K 
does not depend on K . Hence, the sum of all four written out pairs equals to zero. We 
get 

dH_ ( r\l r\m r\l rym\ _i_ ( dH r\l 9H r\l \ 

+h ( & - %) + h mq™ - = o. (52) 

V ax A ox^ J \ at m ^ M ot m ^ A J 

Using the definition of the operator L and grouping in the due order the summands 
we can rewrite the formula fl52l) in the form 



Q l MH) + HL X {Q\) - Q\L,(H) - HL,(Q l x ) = 0. 

Since L is the differential operator of the first order we can use the Leibnitz rule for 
differentiation of product. Hence, the last formula turns into 

L x (QlH) - L„(Q[H) = 0. 
Taking into account the definition Q % a = n ^J a , we get the final formula 

Lx(q% = W)a, (53) 
which is the natural generalization of the potential condition of the vector of momentum. 

The subsequent calculations follow the same pattern. We introduce the results abbre- 
viating a little details of justifications. 

Calculate coefficient of dx^dx x dx p dt I ^ lm ^ at the left hand side of the formula (I46p . 
The first summand -^;dt m standing before the first sum ff^6l) gives the six members: 

-QTQ l ,Q r P + Q P Q l xQ r P - Q™Q l P QD- { ) 

The second summand j^dx^ standing before the first sum f )46|) gives also the six 
members: 



(-i)' {§^(Q l xQ r P - Q l P Ql) + B(Q l P Q» - QW + s? - QW) ■ ( 55 ) 

The summand H^^dx x standing under the sign of the second sum (H6T) gives twelve 
members: 
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H d -&Ql + H ^Q l p ) - { Hd -&Q r P + H^Q 1 ,) + (56) 

(H d -^Q\ + h 9 Sq 1 ,) - (h 9 Sq; + h 9 Sq{) ■ 

QQl 

The summand H-^£dt m standing under the sign of the second sum (H6"|) gives twelve 
members: 

HQ? (0q; + §*q0 - HQ- mQ; + 9 ^q 1 p ) + 

HQ™ l^Ql + 9 ^Q l p ) - HQ™ (§tQ r p + d SQ l x) + (57) 
nn m ( ^kn r -u n l \ — nn m ( n r -i- ^kr> 1 \ 

Gathering together all these summands (I54l)- (I5T1) . and using the definition of the op- 
erator L a , we get 



l,[h(q[q; - q 1 p q\)\ + l x [h(q 1 p q; - q^q;)) + l p [h{q 1 p q\ - q{q^\ = o. (58) 

Let us mark that the expression (1581) does not depends on the multi-index K . That 
is the reason why the right hand side of the formula (158]) is zero. 

To obtain the general formula we consider two multiindices J — {ii, ...i s } C {1, 2, ...n} 
of the order s and H = {a>i, ...a s+ i} C {1,2, ...n} of the order s + 1 . Let us calculate 
the coefficient of the member dx^dt 1 ^ in the form dfl . Let us fix the index \i G H . 

The summand J-^cfa^ standing before the first sum of the formula ( H6l) gives members 
consisted of products of s factors, that was obtained from Q % a dx a where a 6 H\/i. It is 
necessary to choose the index i from the multiindex J , since otherwise would arise the 
factor dt % which does not include into the product of differentials under consideration. 
Each index a£S\/i and i 6 J must be taken once and only once, and moreover, a 
permutation of a pair of indices «j and <x,- leads to the permutation of the differentials 
dxf and dx" , and to the changing of the sign of the coefficient of . Hence, while 
ordering the product of the differentials we get the sign corresponded to the parity of the 
permutation. These demands define exactly the determinant 

detail, (59) 

composed of elements Q l a , where ie J, a 6 H \ fi . To define the sign of this determinant 
it is suffice to find the sign of one concrete summands among its expansion. 

Coefficients of members obtained from the summand j^dt m standing before the first 
sum of the formula (H6"j) give members consisting of products of s factors that is obtained 
from the summands Q™dx a with the upper index m (to avoid the arising of the existing 
differential dt m ). By repeating the previous reasoning we verify that the sum of all such 
differentials has as a coefficient 

v 8H 

E^ det H^vll ; (60) 
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Joining (1551) and (IBDl) . gives 



££=i (S + f^™) det IIQ^II^* AJ = rm x 
E^i^(^)det||g| v ||dx s ^v. ^ 

Let us fix the index \i and gather the summands that are obtained from the members 
^^dx^ of the second sum in the formula (j46j) . Coefficients of these summands are 
adjuncts of the elements Q\ in the matrix HQnwII with the natural cyclic ordering of 
rows and columns. After its multiplication by dQ^/dx^ and after summation we get the 
derivative of the determinant with respect to x^ . One obtains 

E^(det||Qi\JIM^V. (62 ) 

Summands obtained from -^£-dt m standing before the second sum of (jl6|) define the 
derivative of the same determinant over t m . This leads to the following expression 

X><#^ (det HQlgi) dx s dt J \ J . (63) 

The addition of ( )62|) and ( )63|) define the action of the operator L M on the corresponding 
determinant. We unify (1ST]) and (162]) and get 

V 

(L„{H) det WQl^Wdx^dt 1 ^ + (det ||Q| V ||) dx^dt 1 ^) . 

n=i 

These expressions do not depend on the multiindex K , hence, they equal to zero. 
The following theorem was proved 

Theorem 7. The necessary and sufficient conditions for the differential form Q f&Sj) to 
be closed is 

V 

(L»{H) det \\Ql\^\\dx^dt^ J + HL^ (det dx E dt f \ J ) = (64) 

/or any choice S and J . 



10 Connection and Curvature 



To find the connection generated by a field of extremals let us turn to the canonical 
system (|5Tj) . 



n dx 11 _ dH 



kdt™ 9[(g*)']J? 
3i m fef " 
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Consider the variational equations, i.e. the system for the derivatives of solutions with 
respect to a parameter. 



nd[(q 



fl 



Jjm. 



dx } 



Here the parameters f] a corresponds to coordinates of the fiber (for instance, we can take 
the value of the function x(t) at t = to). It is natural to consider if as coordinates 
of the standard fiber of the bundle £ . The variational equations on a solution £(■), <?*(•) 
have the form 



dV fj ir tt tj 



at 
du 



at" 



q*q" 

H XX V + H xq *U 



L q"x 
L xx 



(66) 



Assumption 1. Let U,V be a solution of variational equations IR)6}) for the canonical 
system of the Euler equations, defined on a domain 91 of the space t . 
Suppose that the matrix V is defined and invertible in 91 , 



Consider the matrices W m , m = l...n, which are defined by the formulas W t 
UJ^CVf)' 1 on the domain 91. By differentiating W m we get 



rn 
ftp 



_d_ 



(UV- 1 ) = H xx VV~ l + H^UV- 1 + UV-\H q * x V + H q * q *U)V~ l 



By substituting here the definition of W we get the Riccati equation in partial deriva- 
tives ([T]) for the second variation of the functional ([T]) that corresponds to the Hamiltonian 
H: 



dt m 



H xx + H xq *W + WH a , x + WH * a *W. 



q*x 



(67) 



This equation defines the changing of k -Lagrangian planes with the coordinates 
W™(t) along the extremal surface x(t) (See. [2], [3]). 

Theorem 8. Let C be a smooth foliation on a fibre space £ generated by an invariant 
integral & , and its fibers have diffeomorphic projection on a domain 91 of the space of 
variables t . Let the assumption 1 be fulfilled. 
Then the functions 



' ' pa 



d 

dx p 



K3(m,l) 



dS { dS i dx a 



dti dx a dV 



K 



dS l 

dx^- 



(6* 



defines the solution to the Riccati equation ( ffty - 



Proof. 

The functions 



jjm 



nd(q*y 
kdr] 13 



define derivatives from q* with respect to coordinates of the standard fiber. The matrix 
V gives the operator of differentiation of the mapping from the moving fiber into the 
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standard one of the space £ . In view of assumption 1, this mapping is invertible and its 
invers is defined by the matrix 

(v-y = 

The composition UV^ 1 gives 



d[(q*)XdvP d[(q 



*\t]m 



dr\P dx a dx a 
In view of the formula (J3TT) we get 

BA 

[(O t ] t «(*,^g) = ^ 

Functions A 6 are the different form of the same integrand. Hence 

dA d& 



dqf dq 



The theorem is proved. 
□ 



Corollary 1. In the case (k = 1) i.e. for the Weyl transform [33]) we obtain from (E 
that the Hessian of the action-vector 

w = uv- 1 - db 



dx p dx c 



is the solution to the Riccati equation ( fg7| j 



Let us show that ( 1681) defines the connection on the fibre space £ . With this in mind, 
it is convenient to return from the canonical variables to the Lagrangian ones. Since 
the transform Z is involutori, we have again to make the transform Z of the quadratic 
approximation of the Hamiltonian that generates the variational equation (|66|) 



^((H q * q *U, U) + 2(H q * x U, V) + (H XX V, V)). (69) 
As a new variable we take the derivative of fl69l) with respect to V . 

H q * q *U + H q * x V. 

Then W transforms into 

2)^ = (H q * q *U + Hq^V' 1 = H q * q *W + H q * x . (70) 
As a differential form of connection we consider 

( a = dx a - ^x p dt\ (71) 
The operator of covariant differentiation associated with the form ( 17TT) is 
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V v ,x a = v* w -v^x?. (72) 

The operator of projection of tangent vectors (dt, dx) at a point (t, x) of the bundle 
£ on the fiber has the form 

(dt\ dx a ) H- (0, dx a - fSf^dt). 

Horizontal vectors of the connection V are vectors (dt, dx) that belongs to the kernel 
of that operator: dx — tyixdt 1 = . Hence, the horizontal component of the vector (dt, dx) 
is 

(dt^ixdtf). (73) 
A commutator of matrices A and B we will denote, as is customary, by [A, B] . 

Theorem 9. The tensor of curvature of the connection V equals 

&9% dab r 

Proof. 

Let us recall that the exterior covariant derivative of the form of connection ( is called 
the form of curvature of the given connection [5] . The exterior covariant derivative DC, is 
the value of the exterior derivative d( on the horizontal components of vectors (dt, dx) . 
Let us calculate it 

d( = -ffffdx? A dt m + (J^-V)?pdt l A dtf - -^-^df A dA . 

By substitution of the horizontal component of the vector (dt, dx) we obtain ( !74|) . 
□ 

A connection with the zero curvature is called flat. 

Theorem 10. Let the assumption 1 be fulfilled. Then the connection V that was gener- 
ated by the field of extremals is flat. 

Proof. 

Let us rearrange the formula flTDl) . 



dV 

2) = (H q , q M + H^V- 1 = -q^V- 1 . (75) 



Substitute it in (OJ. We get 



d 2 V t/-1 , dy_ V -ldV_ V -l _ I d 2 V t/-1 , dV\r-ldV\/--l 
dfdP v " r dt 1 v dP v I dPdP v ^ dP v dP v 

□ 

Hence, it was shown that fields of extremals define in the given chart the flat curvature. 
It is because that the existence of a field leads to horizontal integrable distribution of 
planes. Its integral surfaces give the full system of horizontal sections of the bundle £ . 
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11 Examples 



As an example let us consider the standard Hopf bundle S 3 — > S 2 . Here S 3 is realized 
in the space C 2 as a unit sphere: (21,22; Wi,W2), \z\ 2 + \w\ 2 = 1. Fibers of the Hopf 
bundle are defined as big circles {e^z, e^w} passing through each point (z, w) . The 
central projection tt from the center of the sphere onto the tangent plane P at the point 
(1, 0; 0, 0) (at the north pole of S 3 ) turns P into a 3-dimensional projective space MP 3 
with the coordinates that corresponds to the three last coordinates of the plane P . It 
may be considered as a chart 21 on the north hemisphere of S 3 The metric on the space 
RP 3 . induced by the projection tt , is 

ds 2 = (1 + V 2 + Cf + C2T 2 {(dv? + Ki) 2 ) + K2) 2 .) • 

Fibers of the Hopf bundle projects on straight lines (one of two families of rectilinear 
generators of the set of hyperboloid of one sheet — projections of torus on S 3 ). So, 
we have on P the family of geodesic <p , and through each point passes one and only 
one geodesic. Each geodesic gives the absolute minimum of the length among all curves 
lying in the chart 21 . Let us show that there is no action-function that simultaneously 
synchronizes all the extremals of the set <p . Indeed, fibers of the Hopf bundle are obtained 
by the simultaneous rotation about the same angle ip in planes z and w . The point 
(1,77; Cij C2) passes into the point 



A = ((cosy? + 77 sin <p), (sin if — 77 cosy), (Ci cos y + C2 sin y), (Ci sin V ~ C2 cos y)). 

The image of the tt projection of the point A is obtained by normalization — dividing 
by the first coordinate. So, we get 

/ sin y — 77 cos (p £1 cos y + ( 2 sin <p £1 sin <p — ( 2 cos ip \ 
\ cos if + 77 sin if ' cos if + 77 sin ip ' cos cp + rj sin ip ) 

Tangent vectors X to fibers of the Hopf bundle are obtained by differentiation with 
respect to ip and putting ip — . Hence, X = ((1 + z? 2 ), (—77(1 + £ 2 ), (Ci + ^2)) • The 
transversality condition for n — 1 is the orthogonality condition in the metric induced 
on P by the projection tt . Orthogonal planes to vectors X correspond to the space of 
zeroes of the differential form £ = (1 + rf)drj + (— 77(1 + C2MC1 + (Ci + vCzjd^ . However, 
£ A dC, 7^ and the form £ is not integrable. It is the obstacle to design action-function. 
Consequently, *p does not generate a field of extremals. The corresponding manifold is 
not a Lagrangian one. 

The similar situation takes place for the quaternary Hopf bundle S 7 — > S 4 [1\ . It is 
realized as a unit sphere {(z, w), ||-2|| 2 + ||"u;|| 2 = 1} in the 2-dimensional quaternary space 
H 2 . A fiber passing through a point (z,w) is defined as a set of points {az,aw} where 
a runs over the unit quaternary sphere ||er|| = 1 . It is easy to verify that each fiber is the 
central section of the sphere S 7 . Through each point of the sphere passes one and only 
one of these fibers. The projection n from the center of the sphere on the tangent plane 
at the north pole (1,0,0,0;0,0,0,0) of the sphere turns P into 7-dimensional projective 
space MP 7 with coordinates (771, 772, 773; Co> Ci> C2, C3) that corresponds to the seven last 
coordinates of the plane P . The metric on the space MP 7 induced by projection 71 is 
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The projection 7r transfers fibers into the set <p of 3-dimensional planes. Through 
each point passes one and only one plane. Each such plane gives the absolute minimum 
to the functional of 3-dimensional volume in RP 7 in the class of variations lying in the 
considered chart. The normalization relative to the first coordinate defines the parametric 
equation of the planes *J5 : 



= Ho H*> (^M^ 3 ); = Ho ( aw )ji U = 0,1,2,3). 

Note that an element of k -dimensional volume is determined by the length of k- 
dimensional multivector in metric being the tensor k -power of the metric of ambient 
space [7J. The integrand in the situation in question is 



/=(i + r + C 2 



\ 



|/| + |J|=3 



D(a) 



The function / is not a convex one and the Weyl construction is inapplicable. How- 
ever, the function / is convex as a function of Jacobians 



D{a) 

Hence, we can apply the developed theory taking k = 3 . Let us use the transversality 
condition from the section 5. It can be verified that the distribution of planes that are 
transverse to fibers of the quaternary Hopf bundle is not an integrable one. The family <p 
does not give the field of extremals and the corresponding manifold is not a 3-Lagrangian 
one. 

Calculations in this case are much more cumbersome and we do not give it here. 

It is possible to calculate the differential form of connection for the distribution of 
normal planes to fibers of the Hopf bundle and the corresponding form of curvature. It 
may be presumed that we obtain the Chern class of the Hopf bundle. 

It would be interesting to consider exotic 7-dimensional Milnor's spheres for which 
fibers are 3-dimensional spheres that are obtained by the action of unit quaternary a 
on S 7 C H 2 using the formula {az, a h wa^} , where h + j = 1 (see. [22]). Here is an 
additional difficulty connected with the finding of the explicit expression for metric on 
these spheres. 
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